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We study quantum transport in ballistic s±-wave superconductors where coupling between the two bands is 
included, and apply our model to three possible probes for detecting the internal phase shift of such a pairing 
state: tunneling spectroscopy in a N| s±-wave junction, crossed Andreev reflection in a two-lead N| s±-wave|N 
system, and Josephson current in a s-wave|I|s±-wave Josephson junction. Whereas the first two probes are in- 
sensitive to the superconducting phase in the absence of interband coupling, the Josephson effect is intrinsically 
phase-dependent, and is moreover shown to be relatively insensitive to the strength of the interband coupling. 
Focusing on the Josephson current, we find a 0-7r transition as a function of the ratio of effective barrier trans- 
parency for the two bands, as well as a similar phase-shift effect as a function of temperature. An essential 
feature of this s±-wave model is non-sinusoidality of the current-phase relation, and we compute the depen- 
dence of the critical current on an external magnetic field, showing how this feature may be experimentally 
observable for this system. We also comment on the possible experimental detection of the phase shift effects 
in s±-wave superconductors. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.70.Dd 



I. INTRODUCTION 



During the last few years, multiband superconductivity has 
again been at the forefront of condensed matter physics, and 
particularly so after the discovery of high-temperature su- 
perconductivity in the family of intrinsically multiband iron- 
based materials 1 23 . As with all newly discovered supercon- 
ductors with unconventional behavior, one principal question 
is to determine the pairing symmetry of the superconductor. 
In the pnictide superconductors much effort has been devoted 
to this central issue, so far without entirely conclusive an- 
swers. Nevertheless, the leading contender has for some time 
been s±-wave pairing 4 , which in its simplest realization for 
the iron-based superconductors means that the hole-like and 
electron-like Fermi surfaces both host s-wave superconduc- 
tivity, but with opposite sign of the order parameter. (In the 
past, similar sign-shifted order parameters have also been con- 
sidered as a candidate pairing state e.g. of high-T c cuprates^.) 

Distinguishing such a state from an isotropic s-wave pair- 
ing state is highly non-trivial, since both s-wave and s± -waves 
states have the same symmetry, and do not have nodes in the 
order parameter on the Fermi surface. In order to establish 
conclusively the internal phase shift characterizing a possible 
s±-wave state in the iron-based superconductors it is therefore 
crucial to devise phase sensitive pairing probes. A large num- 
ber of proposals for such experiments have been put forth in 
the literature recently. Theories for multiband tunneling spec- 
troscopy have been developed 6 7 8 9 10 11 as well as calcul ations 
of the surface density of states for a s± -superconductoi^ 2 ^. 
In a related context, Andreev bound state s (ABS) a re of- 
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ten pointed out as possible pairing probes 
other class of experiments sug gested inv olves Josephson junc 
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tions, both single j unctiond 1 9 I 20 I 2 1 1 22 ^ , injunction loop d 24 ! 25 -! 
and also various corner geometries employed for Josephson 
interferometryS^El. Yet another work considered possible 
signatures in the AC Josephson effect 28 . In addition, we 



should mention that the Josephson effect for multiband su- 
perconductors with sign-shifted order parameters has previ- 
ously been discussed also in the context of MgB2 29 and bi- 
layer cuprates^l 
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FIG. 1: (Color online) Schematic drawing of the systems under con- 
sideration in this work: a) The model of N|s±-wave junction for 
tunneling spectroscopy as studied in Sec. |III A| b) the model of the 
two-lead N|s±-wave|N junction for the study of crossed Andreev 
reflection in Sec. 
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and c) the model of the s-wave|I|s±-wave 
Josephson junction considered in Sec. |III C] For system (b), we have 
illustrated how an electron in the left-hand lead is converted to a hole 
in the right-hand lead by the formation (together with a electron from 
the right-hand lead) of a Cooper pair in the superconducting inter- 
layer. 

Of the probes listed above, tunneling spectroscopy is prob- 
ably the one that is experimentally most accessible (see Refs. 
|3T1 [32] and references therein), and results here are routinely 
compared with the theory of Blonder, Tinkham and Klap- 
wijk (BTK) for Andreev reflection 33 . Recently, one theoret- 
ical work 9 augmented the BTK-approach to also incorporate 
interband scattering in the superconducting region, which was 
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shown to result in interference effects and subgap bound states 
in the conductance spectra. However, as pointed out soon 
afte^SD, the phenomenological approach employed in Ref. 9 
may fail to capture the effect of interband coupling correctly. 
In this work, we will present an alternative approach of in- 
cluding interband scattering into the BTK framework. 

Another probe which has not been considered in the lit- 
erature so far is crossed Andreev reflection 34 (CAR). This 
is a process contributing to the nonlocal conductance in a 
two-lead normal metal / superconductor junction 35 in which 
an electron impinging on the superconductor from one of 
the leads is converted to a hole in the other lead. This 
phenomenon has previously attracted attention as a possi- 
ble probe both for ferromagnetic superconductors 36 and non- 
centrosymmetric superconductors 3 ^! However, crossed An- 
dreev reflection has not yet, to the best of our knowledge, been 
considered in the context of the s±-wave pairing state. 

Yet another possible experimental signature, which was 
first proposed in the context of iron-based superconductors by 
the present authors in Ref. [38] is 0-7r transition ;] 39 ! 40 !. To 
explain this phenomenon, we draw upon results from Joseph- 
son junctions with ferromagnetic elements. For such systems, 
e.g. a S|F|S junction, the critical current I c switches sign for 
given thicknesses dp of the ferromagnetic interlayer, result- 
ing in non-monotonous dependence of I c on <1f ■ This phe- 
nomenon is ascribed to the junction switching between being 
a (conventional) 0-junction with zero phase difference across 
the junction in the ground state and a 7r -junction, which has 
phase difference ir across the junction in the ground state. Fur- 
thermore, the critical thicknesses dp of S|F|S systems are of- 
ten temperature dependent, which allows for the observation 
of thermally induced 0-ir transitions at T = T 07r as well as 
transitions as a function of interlayer width. 

The possibility of ir -junctions consisting of s±-wave super- 
condu ctors has been mentioned previously in some theoretical 
works 1 2|24|25|4ij m 28 we showed that 0-ir transitions 
were possible in a diffusive s-wave|N|s±-wave junction both 
as a function of temperature and as a function of the ratio of 
interface resistances for each band. The present work is mo- 
tivated by the question of whether these effects persist in the 
ballistic limit, and we perform a complementary, more com- 
prehensive study of the Josephson effects for a simple model 
capturing the essential features of a s±-wave superconductor 
with interband coupling. We find that the 0-tt transition for 
varying ratio of interband resistance persists, but that non- 
sinusoidality of the current-phase relation is significant for 
the present case. For varying temperature we find a some- 
what weaker phase shift effect, which we will relate to the 
more clear-cut 0-7r transition reported for the diffusive case. 
These results for the temperature dependence of the Josephson 
current can be compared with the non-monotonous Josephson 
current between a multigap and a single-gap superconductor 
previously obtained by Agterberg et alW^ 

The outline of this work is as follows. In Sec. |II]we present 
the theoretical framework that is employed to obtain our re- 
sults. This framework will then be applied first t o tunn eling 
spectroscopy of a N|s±-wave structure in Sec. Ill A after 



in a N|s±-wave|N junction in Sec. IIIB The Josephson junc 



tion, to which we will devote the largest share of attention, 



will be treated in Sec. Ill C| The three experimental setups are 
shown schematically in Fig. [T] Some aspects of our model 
and the possible physical realization of the effects found here 
are discussed in Sec. |IV| and we conclude the present work in 
Sec. M 



H. THEORY 

We consider the Bogoliubov-de Gennes (BdG) equations 
for a two-band superconductor with dispersions Sk x mea- 
sured from the Fermi level Ep and gaps Ax, A = {1,2}, 
which read 
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Above, we have used a fermion basis 
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where rjx,k are fermion operators for band A. Considering 
here positive excitation energies E > 0, the solution for the 
wavefunctions ipx is obtained as a generalized BCS expres- 
sion: 



4>x 



where the coherence functions are 




(3) 



u\ = 1 - v\ = v^-IAaI 2 /^), (4) 

while the phases <fix correspond to the broken U(l) gauge 
symmetry of the superconducting state. For the s± state, we 
have 0i — <p2 — 7T. Note that in Eq. ([TJ, no assumptions have 
been made about the pairing mechanism responsible for the 
presence of energy gaps A a in our model, nor of the origin of 
a possible internal phase shift. Our motivation in this work is 
merely to investigate the experimental consequences of such 
a phase shift, when present. 

In order to capture interference effects between the bands, 
it is important to consider carefully the boundary conditions 
in the presence of interband coupling. The above scenario 
corresponds however to a two-band superconductor with no 
explicit coupling between the bands. (Once again, since we 
make no assumptions on the pairing mechanism, the gaps of 
the two bands in Eq. (]T]l might be implicitly coupled through 
two-particle scattering processes, although whether or not this 
would be the case in a microscopic theory will have no conse- 
quences for the present model.) Hopping between the bands 
will be taken into account by adding a single-particle hopping 
term Hh op to the Hamiltonian: 
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which we will turn to the study of crossed Andreev reflection 
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a / dr[r h (r)ril(r) + ri2(r)ril(r)], 
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where rj(r) are fermion field operators in real space, while 
a denotes the hopping parameter. Upon including the stan- 
dard delta-function barrier potential Vo at an interface, one 
may then write down the full BdG-equations in the system. 



Let us, to be definite, consider an N|s± junction, where the 
superconductor occupies the halfspace x > 0. We then have 
H\& = E^, where 
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It is seen that the two bands couple through the interface scat- 
tering as long as a ^ 0, in a simple model which nevertheless 
should be able to capture the main qualitative effects. 

Before turning to applications of this theory, we state the 
resulting boundary conditions for the N\s± junction. For an 
incoming electron from band A' = 1 on the N side (x < 0), 
we write the wavefunction as 

V>jv = [l,0,0,0](e lfex + r 1 e- lfea: )+rf[0,l ! 0,0]e lte 
+ r 2 [0, 0, 1, 0]e~ lkx + r£[0, 0, 0, l]e lkx , (7) 

where k — kp. Here and in what follows, we assume that 
the Fermi level Ep is much larger than (A\, E), such that the 
wavevectors simply read kp — \/2mEp. We also take Ep to 
be the same everywhere in the system, since the effect of any 
Fermi wavevector mismatch (FWVM) can be accounted for 
by adjusting the barrier transparency. Note that although the 
formalism used in Eq. |7) imposes the multiband basis also 
on the normal metal wavefunction, this does not necessarily 
imply that the normal metal has two physically distinct bands. 

For an incoming electron from band A' = 2, the N side 
wave function is simply obtained by letting [1,0, 0, 0]e lfcx go 
to [0, 0, 1, 0]e lkx in Eq. dTJ. Here, {r\ : r^} are the normal 
and Andreev reflection scattering coefficients for band A. We 
let the wavefunction on the superconducting side (x > 0) be 
unspecified for the moment. The general boundary conditions 
can then found from Eq. |6]) as 

i/j n (x = 0) = i) S {x = 0), 
(d x ips ~ d x ip N )\x=o = 2m[Vo diag(i, 1) 

+ a offdiag(l, i)]V»jv, (8) 

where 1 is the 2x2 unit matrix and diag and offdiag de- 
note diagonal and off -diagonal 4x4 block matrices in which 
these unit matrices are embedded. At this point we also intro- 
duce two dimensionless parameters characterizing the system, 
namely the barrier strength Z — 2mVo/k and the interband 
coupling strength a — 2ma/k. 



HI. RESULTS 

A. Conductance spectra 

As a first application of our model, we calculate the con- 
ductance of a N|s±-wave junction and compare it to that of 
its s-wave counterpart. This was also done in Ref. [9l but 
in contrast to their approach, we construct our wavefunctions 
and boundary conditions from the full 4x4 BdG-equations, 
as required for a multiband scenario. In this case, the wave- 
function on the superconducting side reads 

i's = si[u u vxe-^ 1 , 0, 0}e lkx + h[v 1 e 1 ' t ' 1 , Ul , 0, 0}c-' kx 

+ s 2 [0, 0, ua, v 2 e-^]e lkx + i a [Q, 0, v 2 e^ , u 2 ]e' lkx , 

(9) 

with {s\,t\} being the transmission coefficients for band A. 
We will use the gauge <\>\ = 0, and make explicit use of the 
internal phase shift by writing e 1 ^ 1- ^ 2 ^ = S = ±1 for the 
superconductor being a two-band s-wave superconductor or a 
s±-wave superconductor, respectively. For the normal metal 
side, we use ipw from Eq. Q. We can then solve Eqs. ([8]) 
for the given wavefunctions, but as the resulting expressions 
for {r\, r£,s\,t\} do not allow a simple interpretation in our 
case, we give the solution in Appendix A. 

To illustrate the influence of the interband coupling on 
quantum transport, we have plotted in Fig. [2]the probabilities 
of the various reflection processes for an incoming electron 
from band A' = 1 for the case of a transparent interface. For 
decoupled bands, all electrons are Andreev reflected into the 
same band (for subgap energies), and it is shown how this sit- 
uation is altered for a > in a different manner for a s±-wave 
superconductor and a two-band s-wave superconductor. The 
difference between the s-wave and s±-wave case is reduced 
for increasing Z relative to a, and |r 2 | 2 and \r£ | 2 are in gen- 
eral decreased by increasing Z and increased by increasing 
a. Interband scattering also effectively acts to reduce the in- 
terface transparency, although less so for the s± state. Apart 
from these general relations, the dependence of the probabil- 
ities on the coupling a is by no means trivial, and we do not 
attempt to give any further physical interpretation of this pa- 
rameter. 

The conductance for a two-band superconductor normal- 
ized to the normal state conductance Go may, within the BTK 
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FIG. 2: (Color online) Comparison of the probabilities of the re- 
flection processes in a N|s±-wave junction and a two-band N|s- 
wave junction as described in the text. We have chosen zero barrier 
strength, Z = 0, and gap ratio ta = 1.5. We have used a value 
5 = 1 for the interband coupling (for the s-wave and s±-wave case), 
while for the decoupled case we have a — 0. 



formalism, be given as 

G/G 



with G' x = 1 



„A|2 
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2F 
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electron in band A' (see Appendix A), and F 
where the coefficient is evaluated for |Ai| = IA2I = 0. 

In panel (a) of Fig. |3]we have plotted representative results 
for the conductance spectra for different values of the inter- 
band coupling. We have chosen the ratio between the gaps 
somewhat arbitrarily as r& = IA2I/A1 = 1.5, and have in- 
cluded the limiting case of a — 0, which here simply cor- 
responds to the well known BTK result with a double gap 
structure. Furthermore, for values a > Z when Z is small, 
the interband coupling enforces the formation of subgap peaks 
close to the gap edge which are damped and shifted to lower 
energies for decreasing a. This feature becomes more promi- 
nent when 7*a — * 1 (not shown), which makes it observable 
also for larger Z, although also then in a restricted region of 
parameter space. As shown in panel (b) of Fig. [3] no fea- 
tures of this kind appear in the corresponding model without 
an internal phase shift in the superconductor. In the conduc- 
tance spectra of our model, we do not find the very strong 
low-energy conductance peaks reported in Ref. |9l but rather 
features more reminiscent of those of Ref. [10] which may 
be reasonable since their approach was also based on the full 
BdG equations. 
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B. Crossed Andreev Reflection 

One of the most attractive prospects of CAR is as a realiza- 
tion of nonlocally correlated electron states, see e.g. Ref. [42] 
The CAR process is however often masked by the competing 
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FIG. 3: (Color online) Conductance for a N|s±-wave junction (a) and 
a two-band N|s-wave junction (b) for various strengths of interband 
coupling a normalized on its normal state value, where we have set 
Z = 1 and ta = 1-5. 



process of elastic cotunneling (EC), and it is therefore inter- 
esting to search for situations in which CAR dominates. In 
this section, we investigate how the internal phase difference 
of the s±-wave state alters the nonlocal conductance. 

For the left hand side lead (x < 0), we will use the same 
normal region wavefunction (ifa — > t/>l) as in Eq. and for 
the right hand side (x > L) lead we introduce 

tpR = h[l,0, 0, 0}e lkx + tf[Q, 1,0, 0}e- lkx 

+ i 2 [0, 0, 1, 0]e lkx + t£[0, 0, 0, l]e- ikx . (11) 

For the superconducting interlayer (0 < x < L), we now have 
to rewrite the wave function of Eq. |9]into 

4>s = (sie^ x + s a e-^ a )[ui, v u 0, 0] + (s 3 e^ x 
+ S4 e- I9 ^)Kui,0,0] + ( Pl e l ^ x + P2 e~^ x ) 
[0, 0, u 2 , Sv 2 ] + ( P3 e^ x + P4 e-^ x )[0, 0, Sv 2 , u 2 ] 



where we have introduced the wavevectors 



(12) 



(13) 



for electron- and hole-like quasiparticles, respectively. In the 
normal metal regions we can to a good approximation assume 
equal and constant wavevectors k = kp. In our calculations 
we have defined the Fermi energy by the value E-p / A 1 = 1 4 . 

We then apply the boundary conditions of Eq. |8]l to the 
two interfaces at x = and x = L, which results in 16 equa- 
tions in the variables {r\, r£,t\, Si,Pi}, which are solved 
numerically. 
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FIG. 4: (Color online) Nonlocal conductance through a N|s±- 
wave|N junction for a relatively thin superconducting inteiiayer, 
L — 2 ■ 10 fcp , and for ta = 1.5. The upper panels show the 
probability measure for crossed Andreev reflection, while the lower 
panels show for elastic cotunneling. We have used barrier strengths 
Z = (left) and Z — 4 (right), and a number of values for the 
interband coupling a. 



FIG. 5: (Color online) Nonlocal conductance through a N|s±- 
wave|N junction for a relatively thick superconducting inteiiayer, 
L = 8 • 10 fcj , and for ta = 1.5. The upper panels show the 
probability for crossed Andreev reflection, while the lower panels 
show the probability for elastic cotunneling. We have used barrier 
strengths Z = (left) and Z — 4 (right), and a number of values for 
the interband coupling a. 



Since it would have no physical meaning to measure the 
signal for the (virtual) normal metal bands A' separately, we 
choose to consider the average process probabilities 



^ec = ^£(M 2 + M 2 ), 

A' 

PCAR= j£(|tiT+|^| 2 ), 



(14) 

(15) 
(16) 



as the measure of nonlocal conductance, where J2\> again 
denotes summing over incoming electron bands. 

The nonlocal conductance is then proportional to Pec — 
-Pcar, an d we show the result for its separate contributions in 
Figs. [4] and [5] As is expected for the components to the non- 
local conductance, it exhibits oscillations both as a function 
of energy and of the lead separation L, with decaying subgap 
contributions for increasing L. 

It is seen that for high transparency, interband coupling fa- 
cilitates the CAR process with respect to EC, a result which 
may be readily explained, since the coupling acts as an ef- 
fective scattering barrier. Recall that for zero interface re- 
sistance (and no FWVM or spin polarization), the CAR pro- 
cess is completely absent. This result seems to be somewhat 
stronger for a s±-wave superconductor than for a two-band 
s-wave superconductor (not shown), but Pcar is never sig- 
nificantly larger than Pec- All in all, there are only minor 
qualitative differences to be found for the s±-wave state when 
compared to a more conventional s-wave state, and we have 
therefore not included results for the latter here. 



C. Josephson current 

We now turn our attention to the Josephson coupling be- 
tween two superconductors in a S|I|S junction with multiple 
bands. Below, we shall first consider the case where the right 
superconductor is s±-wave while the left superconductor is 
single-band s-wave, with order parameter A s — \A S \ exp (j> s . 
The strategy is to calculate analytically the Andreev bound 
states at the interface, which carry the Josephson current. 
These states are found by using the boundary conditions 
Eq. §8§ for the wavefunctions in each of the superconduct- 
ing regions. However, since we will find that the interesting 
physics stems from allowing different band transmission, we 
let V diag(i, fj -> V = diag(Vi, V x , V 2 , V 2 ). For later ref- 
erence, we also define rz = Z2/Z1 — V2/V1 as the ratio 
between the effective barrier strengths for the two bands; the 
motivation will be discussed in Sec. |IV] Using an alternative 
parameterization to that in Sec. Ill A we write the wave func- 



tion for the left hand side superconductor as 

i> L - Sl [1, e lA , 0, 0}e- lkx + s 2 [e lA , 1, 0, 0}e lkx 

+ s 3 [0,0,l,e lA ]e- lfca; + s4[0,0,e I/3 %l]e lfa , (17) 

while we for the right superconducting region have 

i/) R = ti[l,^-^,0, 0}e lkx + t 2 [e^ 1+lp \ 1, 0, 0}c- lkx 

+ i 3 [0, 0, 1, 5eW 2 -rt]e lkx + t 4 [0, 0, 6e^ +v \ l]^, 

(18) 

where f3 s — arccos (E/\A S \) and (3\ — arccos (E/\A\\). 
The gauge invariant phase difference between the two super- 
conductors has been defined as ip — <\>\ — <f> s . 
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Setting up the boundary conditions of Eq. ([H} yields a sys- 
tem of equations on the form 

At = 0, (19) 

where t = {s\, t\, s!j, t\, s\, t\, sf , t 2 } and A is a 8 x 8 ma- 
trix. The Andreev bound states are found by requiring a non- 
trivial solution for the system, det(A) = 0, which in general 
results in four energy states Er((p). The Josephson current is 
found in the ordinary way by=^ 

* = 2e £f^). (20) 

where Ei denotes the four ABS and f(E) is the Fermi-Dirac 
distribution function. We will define the critical current J c as 
the maximal current allowed by the current-phase relation for 
a given set of parameters, I c = max{I((p)}. We also intro- 
duce the quantity Iq = 2e|Ai| used for normalization of the 
current. 



1. 0-ty phase shifts for varying barrier strengths: the case of equal 
gap magnitudes 



we have for simplicity assumed that \A\\ = |A S | = |A|. To 
show this mechanism explicitly we proceed analytically in the 
limit of a = 0, and this minimal model also serves as a review 
of the basic physics involved in a ballistic Josephson junction. 
Now, the solutions for Eq. (jT9j can be shown to be 

Ef = ±|A| v /l- J D 1 sin 2 (<^/2), (21) 
Ef = ±|A|0 -D 2 cos 2 (<^/2), 

where D\ = 4/(4 + Zf). Ef are the well known solutions 
for a one-band s-wave|I|s-wave junction 43 , while E 2 are the 
corresponding solutions for the negative-gap band. Expand- 
ing to first order in D\ and inserting in Eq. ([20} yields the 
Josephson current 

I = 1 1 sin ip, (22) 

where 1\ — (D\ — D2)Iq/A, It is obvious that for Z2 < Z\ 
one will have D2 > D-y and 1\ < 0, i.e. the system be- 
ing in the ir state. As shown in Fig. |6j the crossover point 
above which the A = 1 contribution dominates instead is 
rz = 1. However, inspection shows that the current does not 
vanish entirely at the crossover point, a fact which is readily 
explained by going to the second order expansion of Eq. ( |21~| >. 
In the limit Z 2 — > Z\ partial cancellation of the two 1st order 
terms then reduces the current to 



0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



s s-t-wave 




I = h sin(2<^), 



(23) 



6 0.7 O.f 



0.9 1 1.1 1.2 1.3 



FIG. 6: (Color online) Critical current for a s-wavejljs±-wave 
Josephson junction as a function of the ratio rz between the effective 
barrier strength of band 2 and band 1 . There is no interband coupling 
(a = 0), and we have set Z\ = 6 and T = 0. Also shown is the 
critical current in the case that the right hand side is a (two-band) 
s-wave superconductor. 

Before investigating the temperature dependence of the crit- 
ical current, we will consider the limit T — > 0. Our main re- 
sults in this section is the observation of a 0-tt transition in the 
Josephson current for varying the barrier strength ratio rz, as 
shown in Fig. [6] For a s-wave|I|s±-wave Josephson junction 
this can be understood in a very simple manner as the compe- 
tition between the A = 1 and the A = 2 band components of 
the current; the band with order parameter Ai = | A. |e 1<?il will 
favor the conventional 0-junction whereas at the same time the 
other band with A2 = — jAje 1 ^ 1 will favor a ir -junction. Here, 



where I 2 = — IqD\/1&. In other words, the second har- 
monic component to the current appears, and is dominating 
close to the transition point. The general non-sinusoidality of 
the current-phase relation close to the transition point is illus- 
trated in Fig. [7] 

Before proceeding, it will be instructive for the subsequent 
discussion to analyze this current-phase relation a little fur- 
ther. In a region close to rz — 1 (and for relatively large 
Z), we may write out the approximate Josephson current to be 
given by the expression 
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sin(2y>). (24) 



For a Josephson junction containing a second harmonic com- 
ponent in the current-phase relation, the ground state needs 
neit her to be a 0-state nor a 7r-state, but may instead be a (p- 



stat d 1 'with a general equilibrium phase difference <po. This 
ground state phase can for our case be found as321 
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(25) 



This phase value evolves smoothly from ipo = ir for < 1 
to ip — for rz ^> 1, passing ip = ir/2 at rz — 1. For 
the case of Z = 6, our model system is a (^-junction for an 
approximate region rz € (0.97, 1.028), and we have verified 
numerically that Eq. ( |24"] > is qualitatively a very good approx- 
imation also well outside this region. The phase difference 
which supports the critical current will on the other hand be 
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denoted as ip*, and can in a similar manner be found to evolve 
from to — 7r/2 for the 7r-state at rz <C 1 to tp* = —ir/4 for 
rz = 1~, where it jumps discontinuously to tp* = 37r/4 for 
rz = 1 + , from which it again evolves smoothly towards ir/2 
for the limiting sinusoidal current-phase relation. This phase- 
shift mechanism will be instrumental to the findings presented 
in Sec. 1111 C 31 
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FIG. 7: (Color online) Current-phase relation for a s-wave|I|s±- 
wave Josephson junction with zero interband coupling (a = 0), with 
Z = 6 and T = 0. 
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FIG. 8: (Color online) Josephson coupling for a s-wave|I|s±-wave 
Josephson junction for various values of interband coupling strength 
a, with a) dispersion of the two Andreev bound states with E < 
shown to the left and b) critical current as a function of barrier 
strength ratio rz shown to the right. Both results are given for in- 
termediate barrier strength Z = 4. 

Next we investigate the general case with nonzero interband 
coupling. The (numerical) solution for the two lower ABS en- 
ergies at zero temperature is shown in panel (a) of Fig. [8]for 
different values of a. While the energy states cross each other 
according to Eq. pT) for a = 0, they repel each other for 
nonzero interband coupling, forming a gap in the ABS energy 
spectrum which increases for increasing a. As a trivial obser- 
vation, this can be understood as a hybridization of the two 
formerly independent bands, as a finite hopping term intro- 
duces off-diagonal matrix elements in A-space. The general 
properties of the current-phase relation remains the same in 
spite of the explicit 7r-periodicity of the ABS dispersion, and 



also here this is explained by the partial cancellation of the two 
ABS contributions. As can be seen from panel (b) of Fig. [8] 
a > does not change the behavior of the Josephson current 
in any dramatic way, and neither does the interband coupling 
influence the position of the Q-tt transition point; it remains 
at rz = 1 for all values of a. This motivates us to suppress 
the interband coupling a in what follows to be able to obtain 
analytically tractable results. 



2. Magnetic field dependence of the critical current 

As a simple application of the model described in the pre- 
ceding section, we now calculate the dependence of the criti- 
cal current I c on an external magnetic field H, i.e. the mag- 
netic diffraction pattern. This quantity is experimentally very 
interesting, and experimental results for I C ( H) have recently 
been presented for iron-based superconductors^M§M2l For our 
model, we are interested in studying how the magnetic diffrac- 
tion patterns depend on the relative barrier strength of the two 
bands, as rz is seen as the primary parameter determining the 
behavior of the system. 

In order to include an external magnetic field to our model 
system, we must define a width W along the z-axis and an ef- 
fective length g?j around x — over which the magnetic field 
H along the ?/-axis penetrates the junction. The magnetic flux 
through the junction is then given by $ = HWdj, and we 
let <l>o denote the magnetic flux quantum. Using the approx- 
imation of Eq. ( f2"4"] > for the current-phase relation, we can 
study our system in the framework of Ref. [45] from which 
we straightforwardly find the expression 
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Evaluating the above expression for the phase difference tp — 
ip* giving the maximum current for the respective rz, we ob- 
tain the Fraunhofer-like diffraction patterns shown in Fig. [9] 
The effect of the 2nd harmonic component to the current is 
evident as a half-integer flux quantum modulation of the crit- 
ical current which grows more pronounced as rz — > 1, but 
whose contribution is vanishing outside the ^-junction region. 
Although results for rz < 1 are not shown here, these are 
largely symmetric with respect to rz = 1. We may also note 
that similar results for the magnetic diffraction were presented 
Ref. 04] albeit for a completely different system. 



3. Temperature dependence of the Josephson current: the case of 
different gap magnitudes 

Motivated by the indications in Ref. 38 that different gap 
magnitudes are necessary for the occurrence of thermally in- 
duced 0-7r-transitions, we now consider a system for the gen 
eral case of A 



III CI 



7^ Ai ^ IA2I. As we showed in Sec 
interband coupling did not affect the 0-7r-transitions as a func 
tion of rz qualitatively, so we will assume in the following 
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FIG. 9: (Color online) Fraunhofer-like magnetic diffraction pattern 
for a s-wave|I|s±-wave junction in an external magnetic field for 
various values of barrier strength ratio. We have used the parameter 
values Z = 6, a = 0. 



that a = 0, an approximation which moreover makes an an- 
alytical approach feasible. Solving the 8x8 system as two 
decoupled 4x4 systems, we obtain the analytical solution as 
given by Eqs. ( |55] l and ( |56| ) in Appendix B. We also refer to 
this appendix for some more information regarding validity, 
existence and uniqueness of this solution. 

We will assume BCS-like temperature dependence of the 
gaps, with the s-wave gap of the left superconductor closing 
at a temperature T CfS = A S (T = 0)/1.76 while both gaps of 
two-band superconductor on the right hand sides close simul- 
taneously at T CjA = T c = Ai(T = 0)/1.76 . We will param- 
eterize the difference in gap magnitudes by r s — A s /Aj and 
Ta = IA2I/A1, and will in most of what follows restrict our- 
selves to r s = 0.5 and r& = 0.3 as a representative set of gap 
ratios, although we stress that our results are valid in a much 
larger portion of parameter space. The resulting temperature 
dependence of the three superconducting gaps is illustrated in 



the inset of Fig. 10 



First, we compare the temperature dependence of the crit- 
ical current both for a s±-wave and a two-gap s-wave super- 
conductor in Fig. 10 The most distinctive feature for both 
these cases is the sharp peak at high temperature. This is 
exactly the temperature T — T* for which two of the gaps 
cross, i.e. |A 2 (T*)| = A S (T*), and although this peak is not 
a signature of the s± -state as such, since it is present irrespec- 
tive of the phase difference between the two right hand side 
gaps, it would be interesting to disclose the mechanism be- 
hind this feature. We turn therefore to the energy dispersion 
of the ABSs, as shown in Fig. [TT]for two temperatures close 
to the peak in the critical current. Firstly, this illustrate how 
E 2 tracks the gap edge of |A 2 (T)| and E\ the gap edge of 
A S (T) < Ai(T) as the temperature is increased, whereas for 
T > T* both states track the smallest of the gaps, i.e. A S (T). 
Secondly, we observe that the energy states are non-dispersive 
for a phase interval centered around ip = and ip = ±7r, for 
Ei and E2, respectively (cf. the discussion in Appendix B), 
so that in these regions the current contributions of the states 
vanish. Thirdly, we also observe that the dispersion of E2 is 



Ai - ■ 


- -|A 2 ] A. 




•<Lo.5 






/' 1 




■< 






/' \ 
jl ] 






1 2 

T/T c , s 


/ ' 1 
J 1 










■ / 

/ 


i\ 
i\ 








s-wave 


i\ 
\\ 
\ 








~s± -wave 


K 


) 


0.2 


0.4 

r/T c , 


0.6 0.8 





FIG. 10: (Color online) Temperature dependence of the critical cur- 
rent for the parameters r s = 0.5, r& = 0.3, with Z — 6 and rz ~ 1. 
Both the results for a two-gap s-wave and a s±-wave superconductor 
are shown. Inset: temperature dependence of gap magnitudes for the 
same parameter set. 



strongly enhanced at T = T* . This last observation can be un- 
derstood by glancing at the expression for E2 in Eq. |56) for a 
given T, from which we realize e.g. by setting cos (<p/2) = 1, 
that the band width of the energy state is at its maximum for 
A s = \A\\. This is of course exactly the case for T = T* . 
Moreover, since the contribution from E\ vanishes for a large 
(^-interval for this temperature, whereas it is non-vanishing 
for E2 for all <p in this limiting case of equal gap magnitudes, 
one does not get the effect of partial cancellation of the two 
current contributions that was present for lower temperatures 
and for A s = Ai = IA2I. We note that although the peak 
strength for these gap ratios is somewhat extreme, we have 
verified that similar peaks or bumps persists in a major part of 
parameter space. 
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FIG. 1 1 : (Color online) Energy of Andreev bound states for two 
temperatures close to the point where |A2(T)j = A S (T) for r 3 = 
0.5 and rA = 0.3. Other parameters are Z = 6 and rz = 1. Shown 
with dotted lines are the relevant gap edges which the energy states 
track. 

We have concluded that a peak in the critical current can- 
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not be taken as a signature of s± -pairing since it results from 
the energy gap crossing of the right hand and left hand super- 
conductor in general. We therefore return to our investigation 
into possible thermally induced 0-ir phase shifts as an unam- 
biguous sign of s±-wave pairing, although apparently no such 
phase shift is present in our results. However, we remember 



from the analysis of the current-phase relation in Sec. IIIC 1 
that in the presence of a second harmonic component to the 
current, a prospective Q-tt transition was smeared out into a te- 
state region for which the critical current remains nonzero. We 
therefore consider the current-phase relation for the junction 
with different gap magnitudes in Fig. 12 for two intermedi- 



ate temperature values. It is evident that the second-harmonic 
component dominates, a fact which can be traced back to the 
vanishing of the ABS contributions for complementary phase 
intervals as discussed above. A related result is that the two 
maxima shown in Fig. 12 originate predominantly from one 
of the energy states each. Furthermore, we have seen that the 
contribution from an ABS is larger the closer the values of the 
gap magnitudes A S (T) and |A A (T)|, and as T -> T* , A S (T) 
and |Aa(T)| are closing in on each other whereas A S (T) and 
Ai(T) are moving apart. Thus the difference in the rate at 
which the gaps decrease causes the E\ state to lose domi- 
nance to E2 for increasing temperature. (Since Ai is by far 
the largest of the gaps, the corresponding ABS dominates for 
T = even though Ai is further from A s than is | A 2 |.) 
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FIG. 12: (Color online) Current-phase relation for a system slightly 
below (T/T c = 0.4) and slightly above (T/T c = 0.7) the thermally 
induced phase shift appearing in the losephson junction with gap 
ratios r s = 0.5, ta = 0.3 and with Z = 6 and rz = 1. The arrows 
indicate the phase difference supporting the critical current (/ > 0) 
for the two temperatures. 



In Fig. 



12 we have also indicated the phase difference ip* 
in the current-phase relation that supports the critical current 
for each of the two temperatures. We now understand that as 
the dominant contribution to the current changes from E\ to 
E2 with increasing temperature, there must be a jump in this 
phase value from <p* > to <p* < 0, and this jump needs to 
happen discontinuously at the temperature T — T v where the 
two contributions balance (cf. our discussion of I c ( r z) in Sec. 
IIIC 1 1. This is our main result in this section: Although the 
Josephson junction is at no point in a 0-state or a 7r-state, the 
system may nevertheless exhibit discernible phase shifts when 
residing in the instate. We illustrate this phenomenon for dif- 



ferent parameters in Fig. 13 and note that similar behavior 
was observed for a large set of different gap ratios as long as 
r s 7^ 1 and r&^l, the basic mechanism behind it being dif- 
ferent temperature dependence of the different gaps. For the 
case of a two-gap s-wave state, a phase shift is of course not 
possible, as the two contributions to the current are then acting 
cooperatively at all times. 
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FIG. 13: (Color online) The upper panel shows the temperature 
dependence of the critical current for the parameters r 3 = 0.5, 
t'a = 0.3 and Z = 6, similarly as in Fig. [10] but for various val- 
ues of barrier strength ratios rz- The lower panel shows the phase 
difference ip* supporting the critical current as a function of temper- 
ature for the same parameter set as above, illustrating the effect of 
the (discontinuous) phase shift from <p* > to p* < 0. 



IV. DISCUSSION 

Comparing the three systems considered in the previous 
section, it is easy to see that role played by interband scatter- 
ing differs fundamentally. On the one hand, tunneling spec- 
troscopy and nonlocal conductance in the absence of inter- 
band coupling is not dependent upon the relative phase differ- 
ence of the two s±-wave order parameters, being merely the 
sum the contribution from two decoupled s-wave states. On 
the other hand, phase information enters explicitly into the 
calculation of the Josephson current, so that the interplay be- 
tween the phases of the two order parameters is evident also 
for zero interband coupling. Furthermore, it seems that the 
behavior observed for the Josephson current remains qualita- 
tively unaltered also for finite a. This explains how it seems 
much more appealing to obtain phase information from multi- 
band superconductors by the use of Josephson junctions than 
by tunneling spectroscopy, and why we will focus our discus- 
sion on this experimental probe. 

To be able to compare our results for the ballistic limit with 
our previously obtained results for the diffusive limit in Ref. 
l38l we now briefly recapitulate this work. Here we employed 
the quasiclassical Usadel equation 50 to study Josephson cou- 
pling in a s-wave|N|s-|--wave junction in the limit of weak 
proximity effect, an approximation which is warranted for 
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low-transparency interfaces. We showed that for this case, 0-tt 
transitions were observed both as a function of barrier trans- 
parency ratio (for arbitrary gap ratios r s , r&) and as a func- 
tion of temperature (for some values of the gap ratios). Here, 
the obtained current-phase relation was purely sinusoidal irre- 
spective of parameter values, a result which can be explained 
by the fact that the linear Usadel equation corresponds to only 
a 1st order approximation in the interface resistance, so that 
no 2nd harmonic terms will appear. Our present results, on 
the other hand, are valid for arbitrary interface resistance, and 
we see that in this model the 2nd harmonic term plays a cru- 
cial role in the behavior of the Josephson junction, which we 
will discuss more below. We should also remark here on the 
difference between the diffusive and the ballistic model in that 
the former in contrast to the latter has a interlayer with finite 
thickness, which was needed to justify the assumption of weak 
proximity effect. 

The importance of a prospective 2nd harmonic contribution 
to the Josephson current is natural when we are concerned 
with 0-7r transitions, as this component may dominate when 
the 1st harmonic component vanishes close to the transition 
point. This fact, and the (^-junction behavior that follows, 
has bee n pointe d out several times in the context of S|F|S 
j unction j 43 l 51 l 52 [ For our model, the influence of the 2nd har- 
monic is seen to be particularly prevalent in the case of differ- 
ent gap magnitudes and/or high interface transparency. Before 
discussing its implication on the thermal phase shift effect ob- 
served here, we consider the case of more conventional 0-n 
transitions. Most often when establishing 0-ir transitions in a 
Josephson junction, one looks for a sharp cusp in the critical 
current as a function of the parameter in question, although 
this method cannot discern which side of the transition repre- 
sents the 0-state and which represents the 7r-state. By using a 
rf SQUID configuration 43 one may however measure the jump 
in the critical phase difference across the junction, which for a 
sinusoidal current-phase relation would be from ip* = w/2 to 
ip* = — tt/2 or vice versa. (Note that it is crucial to this argu- 
ment that one considers a current-biased experiment in which 
a current / > is forced through the junction, letting the 
phase difference adjust accordingly.) In the presence of higher 
harmonics in the current-phase relation the jump from ip* > 
to ip* < or vice versa will in general be different, cf. the 
transition for varying rz, but the principle remains the same. 
This also holds when the sinusoidal component to the current- 
phase relation is subdominant for all parameter values, such 
as for the thermal transitions reported here, so that the critical 
current is not even close to zero at the transition point. In- 



specting Fig. 13 one sees that the critical current does in fact 



reach a minimum at T = T v . Hence this phase-shift effect 
can be regarded as a degenerate form of 0-7r transition which 
can only be established by SQUID measurement of the critical 
phase difference. Alternatively, one could of course demon- 
strate the transition by using SQUID to map out the entire 
current-phase relation 53 , but observing a single phase-shift of 
the critical phase tp* may be simpler experimentally. 

Considering then the peak phenomenon described for the 
temperature dependence of the critical current, as pointed out 
earlier, it does not pertain to the s± -state per se, but is a gen- 



eral result in this framework of two gaps crossing at a cer- 
tain temperature. In fact, this even holds when none of the 
two superconductors are multiband superconductors. Exper- 
imentally, this can however be understood to be a somewhat 
artificial situation, as the phenomenon would not occur for a 
junction consisting of two conventional superconductors with 
different zero-temperature gap magnitudes because of the uni- 
versal ratio 2A(0) /T c = 1.76 for BCS superconductors. And 
for e.g. high-T c cuprates, for which the corresponding ratio 
is larger, the value of A(0) is typically much larger than for 
any conventional superconductor as well. So although multi- 
band superconductors are not necessary as such, the described 
situation can occur here much more easily because the su- 
perconducting pairing for both bands typically vanish at the 
same critical temperature, whereas the gap ratio r/\ ^ 1. This 
is the situation for the conventional multiband superconduc- 
tor MgBjISl, and also seems to be the case for the iron-based 
superconductors 32 . We should note that similar behavior was 
not found in the diffusive case 38 , but that a finite temperature 
maximum in the critical current for multiband superconduc- 
tors was predicted in Ref. 29 In that case, the effect was how- 
ever ascribed to thermal effects combined with different sign 
of the two order parameters, and is not related to gap crossing 
irrespective of the order parameter sign as in our case. Fur- 
thermore, in Ref. |29]as well as in our results for the diffusive 
case, the current-phase relation was implicitly assumed to be 
purely sinusoidal, which may explain some of the differences 
with our present results for the ballistic case. 




FIG. 14: (Color online) Illustration of the physical interpretation of 
the barrier strength ratio rz- The outer circle represents the Fermi 
surface of band 1 of the s±-wave superconductor and the inner circle 
represents band 2. The two dashed circles in between represent the 
Fermi surface of the s-wave superconductor on the left hand side of 
the junction for two situations: largest FWVM with respect to band 
1 (Z\ > Z 2 ) and largest FWVM with respect to band 2 (Z\ < Z%). 



As the dependence of various observable quantities on the 
barrier strength ratio rz was considered frequently through- 



out Sec. Ill C we would now like to present a more thorough 
rationale for this parameterization. Firstly, we note that al- 
though our model assumes the same Fermi wavevector kp for 
all bands in all regions of our setups, any FWVM between the 
different regions is equivalent with an increase in the barrier 
strength Z. And for different Fermi wavevectors fep, s , k-p : \ 
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for the s-wave superconductor on the left hand side and band 
1,2 of the s±-wave superconductor on the right hand side in 
the Josephson junction, respectively, this gives rise to differ- 
ent effective barrier strengths Z\^, for the two bands. This 
line of thought is illustrated in Fig. 14 The idea of tailoring 



the experimental setup by the use of materials with appropri- 
ate Fermi surfaces was first proposed in Ref. [19] and as we 
discussed in Ref. [38] it might be possible to produce a series 
of junction samples with different barrier strength ratios rz 
by varying the doping level in the non-s±-wave region of the 
junction. In this manner, it is conceivable that a 0-ir transition 
can be observed for varying doping level, analogously as to 
how 0-7r transitions are observed S|F|S junctions for varying 
interlayer thickness. The preceding argument naturally raises 
the question whether it might be more appropriate with a pa- 
rameterization in which increased transmittance through the 
barrier for one of the bands was accompanied by decreased 
transmittance for the other band. We did nevertheless define 
rz as simply the effective barrier strength ratio because it is 
hard to tell exactly how the relative transmittance will change 
with doping level. This does probably also make it quite chal- 
lenging to experimentally produce the right series of samples 
to observe 0-7r-transitions, which is what makes possible ob- 
servable signatures for varying temperature all the more ap- 
pealing. 

Finally, we discuss our model in context of the recently 
discovered iron-based superconductors. It should be stressed 
that our model is to be taken as a minimal model describ- 
ing the generic behavior of transport phenomena in a two- 
band s± -superconductor, but we would like to point out how a 
more realistic model should take into account the specifics of 
the iron-based superconductors. Firstly, the BTK approach 
does not incorporate any details of the band structure, and 
spherical Fermi surfaces are assumed. Secondly, ours is a 
two-band model whereas it has been argued that at least four 
bands, two hole-like (h) and two electron-like (e), should be 
included to capture the physics behind the superconductivity 
in these materials 55 . The main effect of their inclusion from 
the point of view of transport would be the possibility of e-e 
and h-h interband scattering between nearly-degenerate elec- 
tron and hole bands, respectively. Since these scattering pro- 
cesses would involve no internal phase-shift, we expect that 
the result would be qualitatively similar to the two-band case. 
This assumption can also be justified by the fact that h-h and 
e-e scattering processes should be weak in iron-based super- 
conductors compared to the spin density-wave-enhanced e-h 
interband processes^, so that we to a good approximation can 
consider degenerate e and h bands. Generalization to a non- 
degenerate four-band model could nevertheless be made in 
our theory by a straightforward extension of Eqs. ([l} and Q, 
with the inclusion of h-h and e-e interface scattering terms in 
Eq. |5]), although an analytical treatment in that case would 
be a daunting task. Furthermore, one might have gap magni- 
tudes that were momentum dependent, but the approximation 
of constant s-wave gaps on each of the Fermi surfaces should 
be reasonable. (The possibility of a d-wave gap or other pair- 
ing symmetries with nodes on the Fermi surface is left out of 
the question in this work, since the majority of experiments 



so far seem to indicate a nodeless gap on the Fermi surface.) 
Another extension would be to include interband scattering in 
the bulk of the s± -superconductor, and not only near the in- 
terfaces as in our case, or even more sophisticated models e.g. 
including momentum dependence in a. 

Regarding the magnetic field dependence of the critical 



Josephson current described in Sec. IIIC2 we may compare 
our results with the experimental results for iron-based super- 
conductors available at the moment. Inspecting the diffraction 
pattern in Fig. 3 of Ref. [48] we note an intriguing similarity 
with ours for rz > 1 in that the critical current is nonvan- 
ishing between the diffraction maxima. This may however 
just as well be the combined result of nonuniform current dis- 
tribution, trapped flux and deviation from the small junction 
limit 57 , so that we cannot with any certainty interpret this ob- 
servation as evidence for a non-sinusoidal current-phase rela- 
tion, nor would non-sinusoidality necessarily imply s±-wave 
pairing. (The diffraction patterns obtained in Refs. [47] and 
|49|can on the other hand not be compared with our results at 
all, as the experimental situations for those works are differ- 
ent.) It should also be noted that our modelling of the flux 
threading the junction is rather simplified, and does not in- 
clude effects that may be present in real samples 57 . More 
importantly, assuming isotropic order parameters and Fermi 
surfaces, our model is insensitive to how the junction geom- 
etry is chosen. We therefore cannot capture the directionality 
of the electron-like Fermi surfaces in the folded Brillouin zone 
of iron-based superconductors, which is essential in other pro- 
posals for phase-sensitive corner junctionals anc j related ge- 
ometries. 

It would also be very interesting to see how robust the 
results presented here are to the introduction of material 
impurities. The iron-based superconductors are mostly ex- 
pected to reside in some intermediate regime of impurity 
concentration 2 , thereby making neither the ballistic nor the 
diffusive limit a completel y accurate de scription. In fact, a 
number of theoretical works^ 5 ^ 6 ! 17 ! 18 ! 58 ! 5 ^ depend upon a sig- 
nificant influence of impurities to explain the experimental re- 
sults or to induce experimentally observable bound states. Our 
study in Ref. [38] was motivated by the fact that the diffusive 
regime is often the experimentally relevant one. Although tak- 
ing the diffusive limit may not be strictly valid in this case, the 
results found might nevertheless capture important features of 
the real materials. In light of this, it would be very interesting 
to compare the results obtained in the diffusive and the ballis- 
tic limit with calculations performed using the quasiclassical 
Eilenberger equation 60 , which allows for arbitrary impurity 
concentration. This would require a multiband extension of 
the Zaitsev boundary conditions 61 , and such a theory has only 
very recently been developed (see Ref. |62| |. 



V. CONCLUSION 

Possible signatures of s±-wave pairing in tunneling spec- 
troscopy stem mainly from the multigap nature of the super- 
conductor, but also from interference effects when the inter- 
band coupling is strong relative to the barrier strength. This 
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may lead to subgap peaks in the conductance spectra not 
present for a corresponding s-wave model, although the ap- 
pearance of these are relatively sensitive to the parameter val- 
ues used. Similarly for the nonlocal conductance, it is found 
to be very difficult to discriminate qualitatively the interfer- 
ence effects of a s±-wave state from those of a two-band s- 
wave state. Josephson coupling is on the other hand an in- 
trinsically phase-dependent phenomenon, so it is natural that 
it is here that we find the most promising signatures of s±- 
wave pairing, namely Q-tt transitions or similar phase shifts 
in a s-wave|I|s±-wave junction. These are neither dependent 
on, nor considerably affected by, the presence of interband 
coupling. As in the diffusive case®, we find 0-n phase shifts 
as a function of the relative interface transparency, an effect 
whose detection is possible in principle, but difficult in prac- 
tice. We have also shown that a phase-shift effect is present 
as a function of temperature, and although this effect is not as 
robust as the one reported for the diffusive case, it may never- 
theless be possible to observe using a SQUID setup. For both 
cases, we have shown how the phase shifts can be ascribed 
to the competition between Andreev bound states for the two 
bands, and how the non-sinusoidality of the Josephson current 
is essential in the description of the phase shifts. We have also 
pointed out that this 2nd harmonic component in the current- 
phase relation may induce half-integer quantum flux modu- 
lations in the magnetic diffraction pattern of the Josephson 
junction. In addition, we found a peak feature in the tempera- 
ture dependence of the critical current for the case of different 
gap magnitudes, an effect ascribed to the crossing of two gaps. 
Although it is hard to tell how relevant the signatures reported 
in this simplified model are for possible experimental realiza- 
tions of the s±-wave pairing state, our results shed more light 
on the basic mechanisms of transport and their implications in 



transmission coefficients given by 
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For the case of A' = 2, the corresponding expressions read 
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The reflection coefficients are then found for both cases by 
insertion into 



ri = -<5 A ',i + U1S1 + vih, 
rf = V1S1 + U\t\, 
T2 = —Sx',2 + u 2 s 2 + 5v 2 t 2 , 
r% = Sv 2 s 2 + u 2 t 2 , 



(35) 
(36) 
(37) 
(38) 



where 8y_i is the Kronecker delta. 

The auxiliary quantities used for A' = 1, 2 are given by 
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APPENDIX A: REFLECTION AND TRANSMISSION 
COEFFICIENTS FOR THE N|s±-WAVE JUNCTION 



In this section, we give the analytical solution for the re- 
flection and transmission coefficients. We have to consider 
the two cases A' = 1,2 for the incoming electron band inde- 
pendently, but will use the same symbols for the coefficients 
to simplify notation. First considering A' = 1, we have the 



A = u%U2 - 5viv 2l 
B = viu 2 - 6U1V2, 
C\ = v\ - u\, 
7A = 4u 2 - C X Z 2 . 



(50) 
(51) 
(52) 
(53) 



The expressions above are valid both the s±-wave and the 
coupled s-wave case, where s±-wave is found by setting 
S = — 1 and s-wave by S = 1. 
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APPENDIX B: SOLUTION FOR THE ABS ENERGIES FOR the general case of different gap magnitudes yield after some 
DIFFERENT GAP MAGNITUDES manipulation the equation 

The coefficient matrix A for each of the uncoupled bands in 



9m|A A | = (4 + Z\) sin (20„ + 2/? A ) - Z\ [sin 2/3 s + sin 2/3 x ] - 85 x sin (ft + X ) cos ip = 0, 



(54) 



with Si = 1 and 5 2 = —1 for a s±-wave superconductor. Using that cos j3 x — E/\A\\ and cos f3 s — |A A |/A S cos /3 A we can 
solve the equation for cos 2 j3\, which produces the solutions 



Ai sin 2 (ip/2) + Bi- JCi sin 4 (ip/2) + D l sin 2 {ip/2) + F l 
2{Z\ + 4) 



. , Z\ + 2 / A 2 cos 2 {ip/2) + B 2 - ^C 2 cos 4 {ip/2) + D 2 cos 2 (p/2) + F 2 

Ejc, =± \ 



2(Z 2 +4) 



(55) 
(56) 



in addition to several other unphysical solutions. The auxil- 
iary quantities here are given by 





Ax = 2Kx, 


(57) 


B x 


= A 2 s + \Ax\ 2 -K x , 


(58) 




C x = 8#aA s |A a |, 


(59) 


D x 


= 4{A s -\A x \) 2 Kx, 


(60) 




F A = (A S -A 2 ) 2 , 


(61) 


K x = 


8A s |A A |/(Z 2 + 2) 2 . 


(62) 



To justify that the given solutions are the only solutions and 
are also in fact valid for all parameters, we have verified nu- 
merically that Rc|A A | = 3m|A A | = for all solutions of 
E^ used in this work. However, as can be seen by compar- 



Ex {(p) for around ip w for Eq. ( |55| ) or around <p w ±7r for 
Eq. ( |56| ) does not produce a valid result for |A A | ^ A s . The 
explanation is that the physical Andreev bound states simply 
vanish in these regions, and we have again confirmed numer- 
ically that |A A | = have no real solution for E here. In fact, 
solving only for the imaginary part of the determinant yields 
(clearly unphysical) solutions \E\\ > min{A s , |A A |}, which 
furthermore result in complex factors sin/3 A , rendering Eq. 
p4] ) invalid as an expression for the imaginary part of the de- 
terminant. In the results presented above, we have handled 
this numerically by setting the bound state energy equal to 
the gap value when vanishing, so that it does not contribute 
to the current (since the energy states vanish at the gap edge 
with zero slope), although the energy states strictly speaking 
do not exist at all in these regions. 



ing with Fig. 1 1 and the accompanying discussion, evaluating 
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